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Abstract 

We demonstrate that the conventional application of linear models to the analysis 
of optoelectromechanical properties of nanostructures in bandstructure engineering 
could be inadequate. Such linear models are usually derived from the traditional 
bottom-up approach applied to the analysis of nanostructure properties. At the 
same time, in the hierarchy of mathematical models for semiconductor device mod- 
elling constructed on the basis of the top-down approach, we deal predominantly 
with models where nonlinearity is essential. In this contribution, we analyze these 
two fundamental approaches in bridging the scales in mathematical models for the 
description of optoelectromechanical properties of nanostructures. The focus of the 
present paper is on a model based on the coupled Schrodinger-Poisson system where 
we account consistently for the piezoelectric effect and analyze the influence of dif- 
ferent nonlinear terms in strain components. The examples given in this paper show 
that the piezoelectric effect contributions are essential and have to be accounted for 
with fully coupled models. While in structural applications of piezoelectric mate- 
rials at larger scales, the minimization of the full electromechanical energy is now 
a routine in many engineering applications, in bandstructure engineering conven- 
tional approaches are still based on linear models with minimization of uncoupled, 
purely elastic energy functionals with respect to displacements. Generalizations of 
the existing models for bandstructure calculations are presented in this paper in the 
context of coupled effects. 
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1 Introduction 



In low-dimensional semiconductor nanostructures (LDSN) the motion of elec- 
trons can be confined spatially, from one, two, and even three spatial direc- 
tions. In the latter case, such nanostructures are known as quantum dots 
and often termed by physicists as O-dimensional structures, reflecting the fact 
that the motion of carriers is constrained from all three spatial directions. 
These structures have been receiving an increasing interest due to new tech- 
nological advances and fascinating applications they offer. Indeed, they can be 
used as biological tags in cell biology and biomedicine, be used in construct- 
ing quantum bits for quantum computing, and be applied in a wide range of 
more traditional structure- and device-like applications, including photodec- 
tors, laser-based emitters, etc. 

While in many such applications the focus is on optical properties of future 
devices, it is important to remember that the formation of LDSNs, and in 
particular quantum dots, is a competition between the surface energy in the 
structure and strain energy. Hence, mechanical properties are essential in de- 
signing quantum-dot-based devices and structures. Further, many quantum 
dot structures have a well pronounced piezoelectric effect which does con- 
tribute to their overall properties in a non-trivial manner. These coupled elec- 
tromechanical effects will become increasingly important for the current and 
future applications of such nanostructures. In designing stable strained nanos- 
tructures, computational modelling provides a major tool for predicting their 
optoelectromechanical properties. 

During the last decade, the attention of the science and engineering community 
to the influence of strain effects on quantum mechanical properties of LDSNs 
has been growing rapidly [1,2,3,4,5]. A majority of the published works were 
focusing on strain effects only, without taking into account electromechanical 
interactions due to the piezoelectric effect. Those authors who did account for 
the piezoeffect based their considerations on the minimization of uncoupled, 
purely elastic energy functionals with respect to displacements. Under this 
approach, the Maxwell equation for piezoelectric solids and the equations of 
elasticity were effectively solved in either uncoupled or semicoupled manner. 
Pan with his collaborators [6,7] were the first who have attracted the attention 
of the bandstructure engineering community to the importance of coupled ef- 
fects. Based on his semi- analytical Green's function approach applied to an 
idealized half-space structure [6], he demonstrated that only the fully cou- 
pled model can lend a reliable prediction. Further, based on a combination of 
analytical (for the ID case) and numerical (for the 2D case) techniques, the 
idea was generalized to the device level [7], but no details of the numerical 
procedure were given. All the above studies were based on the linear theory 
only. Furthermore, the nature of the analyzed problem allowed a number of 
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simplifications, including those related to the wetting layer. The inclusion of 
the wetting layer in a consistent manner not only increase the computational 
complexity of the problem in several times due to different spatial scales, but 
may also require the formulation of non-trivial boundary conditions [8] . 

In this paper, we base our consideration on the coupled Schrodinger-Poisson 
model where we account consistently for the piezoelectric effect and analyze 
the influence of different nonlinear terms in strain components. We structure 
the paper as follows. In Section 2, we analyze two fundamental approaches 
in bridging the scales in mathematical models for the description of optoelec- 
tromechanical properties of nanostructures. In Section 3, we provide the core 
model for the description of coupled electromechanical interactions in piezo- 
electric semiconductor solids. The model is exemplified for hexagonal (WZ) 
and cubic (ZB) materials used in our computational experiments. In Section 
4, we give details of the model for bandstructure calculations, focusing on the 
k • p approximation as a convenient framework for incorporating strain and 
piezoelectric effects. A general procedure for modelling quantum dot nanos- 
tructures, based on the variational formulation of the problem, is outlined in 
Section 5. In Section 6 we provide details of numerical experiments demon- 
strating the influence of the piezoelectric effect and analyzing contributions of 
nonlinear terms in strain components. Conclusions are given in Section 7. 



2 Waves propagation in anisotropic media: applying experience 
from solid and fluid mechanics to bandstructure engineering 



Waves propagation in anisotropic media has always been a topic in the heart 
of scientific inquiries and a source of new ideas for engineers. This topic is of 
immense practical importance in the context of both solid and fluid mechanics. 
Purely elastic waves in shells and other structures have been studied in the 
context of fluid-solid interactions at least since the late 1950ies, providing 
many meaningful examples where coupling effects become essential. Around 
the same time, an increasing interest to coupled problems was also generated 
by the analysis of dynamic thermal stresses in structures, in particular by the 
famous Danilovskaya problem, first formulated in 1950. Since then, coupled 
effects, in particular in anisotropic materials, have continued to fuel interest 
to their studies due to both, theoretical challenges and an increasing range of 
practical applications. 

In what follows, we will focus on the coupling between electric and mechan- 
ical fields in low-dimensional semiconductor nanostructures. The core of the 
mathematical models dealing with this sort of coupling contains the equations 
for coupled electromechanical motion of piezoelectric solids. The experience 
accumulated in mechanics of solids in solving such equations becomes now 
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invaluable in the area known as bandstructure engineering and in the mod- 
elling semiconductor quantum structures in general. At the same time, the 
experience accumulated in fluid dynamics applications is equally important 
in this area. Indeed, the hydrodynamic approach in the analysis of semicon- 
ductor devices, superlattices, and other semiconductor structures has been 
an important tool in semiconductor modelling for a number of years. Fluid- 
dynamics-like hydrodynamic approximations have been widely utilized in this 
area and more recently several their extensions have been proposed to account 
for quantum effects. The importance of such semiconductor systems as quan- 
tum wells, wires, dots, and superlattices [9] will continue to grow in nanoscale 
electronics, photonics, and bioengineering. New technological advances in ap- 
plications of these structures require to have a fresh look at their modelling 
aspects, in particular in the context of their optoelectromechanical properties. 
While strain effects are fundamental to such properties, in the bandstructure 
engineering literature their influence is still typically analyzed with simplified 
linear models based on the minimization of uncoupled, purely elastic energy 
functionals with respect to displacements. The applicability of such models 
is limited as coupled effects related, e.g., to built-in spontaneous and piezo- 
electric polarization become essential. New models accounting for these effects 
need to be developed. 

The modelling experience accumulated in both mechanics of solids and fluid 
mechanics can help in achieving this task. To get started, note that in semi- 
conductor systems we are dealing with, both classical and quantum effects 
are interlinked, and the analysis of such systems and the choice of modelling 
tools depend critically on the spatio-temporal scales required for specific ap- 
plications. We can start constructing a model for the analysis of such systems 
from the fundamental quantum level by specifying the Hamiltonian of the 
system within the Schrodinger framework. However, then we should incorpo- 
rate additional effects, pronounced at larger scales, such as piezoelectric, into 
the obtained approximate model. This is the bottom-up approach to modelling 
semiconductors, applied actively today for the analysis of nanostructures. Al- 
ternatively, we can attempt to carry out some physics-based averaging right 
from the beginning, applying the top-down approach to modelling semiconduc- 
tors. A good example, clarifying conceptually the applicability and limitations 
of these approaches, can be provided by considering models for superlattices. 
Based on the underlying physical assumptions, there are two major classes 
of such structures, classical and quantum. These structures have additional 
periodicity on a scale larger than atomic. The idea of creating quantum su- 
perlattices is due to L. Keldysh (1962). Experimentalists reported the creation 
of such objects about a decade later (L. Esaki, 1970; Zh. Alferov et al, 1971). 
It has been the domain of solid state physics where tools of solid mechanics are 
essential and well established. On the other hand, the idea of creating classical 
semiconductor superlattices was originated from a fluid mechanics analogy. It 
is well known (e.g., [10] and references therein) that a fluid under gravity can 
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reach its equilibrium if its temperature depends on the height only. If the tem- 
perature gradient, directed down, exceeds certain critical value, we observe a 
free convection of the fluid. If we assume that this process takes place between 
two infinitely long horizontal planes heated to different temperatures (tem- 
perature of the lower plane is higher), under a sufficiently high temperature 
gradient the fluid becomes unstable and we observe stationary convective mo- 
tion. Due to the underlying assumptions, in the horizontal plane the motion is 
expected to be periodic. Based on this hydrodynamic analogy, a similar idea 
was proposed in early 1970ies in the context of carrier motion in semicon- 
ductors, where the role of gravitational field could be played by the electric 
field with heating produced by, e.g., light ([10] and references therein). The 
analogy between wave phenomena in classical superlattices and the behaviour 
of wave functions of an electron moving in a periodic potential field of a quan- 
tum superlattice can be exploited when developing a hierarchy of models for 
the analysis of semiconductor structures. In both cases, we have an additional 
periodicity of the structure. However, the difference between these two cases 
lies in the fact that while for classical superlattices such an additional peri- 
odicity leads to the quantization of wave energy, in quantum superlattices it 
leads to the quantization of carrier energy. If the period of the potential field 
exceeds the length of free carrier runs, so that on this specific spatio-temporal 
scale carriers will not be affected by the action of the additional periodic field, 
semiconductor superlattices behave like classical structures. In the latter case, 
many (fluid-mechanics analogy based) techniques developed for semiconduc- 
tor device modelling at sub-micron scales can often be applied. The study of 
quantum superlattices requires more fundamental approaches to account for 
atomic scales. 

Note that already at the classical level we have to construct a multiscale hier- 
archy of the models. Indeed, the standard drift-diffusion approximation may 
not be an appropriate modelling tool even for classical superlattices, while 
hydrodynamic and kinetic models provide quite useful tools for the analysis 
of such structures. Based on relaxation time approximations, a classification 
of the hierarchy of mathematical models for these structures was discussed 
in [11,21]. At the top of this hierarchy is the Liouville equation framework 
which leads to substantial difficulties in practical realization of this approach. 
Hence, most practical approaches stem from the kinetic-type models such as 
the semi-classical Boltzmann equation. In these models, scattering of carriers 
on each other is not essential. However, the scattering of carriers on imper- 
fections of the lattice plays the dominant role. Hence, charge carriers under 
this approach cannot be considered as an independent thermodynamical sys- 
tem. Based on the moment methodology or the Hilbert expansion method, 
a range of macroscopic models can be derived, among which hydrodynamic- 
type models play a prominent role. In such models the electron-hole "plasma" 
can be considered as an almost independent thermodynamical system that 
only weakly interacts with the crystal lattice. This group of models as well 
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as quasi-hydrodynamic models can account for non-equilibrium and non-local 
behaviour of semiconductor carrier "plasma". Further details of the devel- 
oped computational techniques for such models can be found in [21,11,12]. 
We note that the equations we deal with in such situations are similar to, 
but differ from, the hydrodynamic equations of fluid mechanics. They consist 
of the Poisson equation, equations of continuity (for carrier concentrations) 
and energy transfer. Furthermore, re-distributions of charge carriers lead to 
an additional field, a phenomenon absent in the fluid mechanics. Attempts to 
apply these types of models to other semiconductor structures and devices at 
smaller scales have led recently to the development of extended hydrodynamic 
models that should incorporate quantum corrections [13]. The important ob- 
servation is that all the models we have discussed above within the top-down 
approach are intrinsically nonlinear. 

At the other end of the spectrum of model hierarchy are the models developed 
with the bottom-up approach. Surprisingly, up to date the majority of research 
efforts in this area has been concentrated on linear models. In what follows, we 
focus on the analysis of quantum dot structures and show that the conventional 
approaches to the analysis of these structures based on linear models need to 
be augmented to account for coupled nonlinear effects. 



3 Coupled electromechanical interactions in quantum dot nanos- 
tructures 



Mechanical effects profoundly influence electronic and optical properties of the 
nanostructures. Two points should be mentioned in this context. Firstly, we 
note that the key to intrinsic properties of quantum dot structures lies with 
strain effects arising from lattice mismatch. Following [14], where the authors 
started their reasoning from the total Helmholtz free energy function, we as- 
sume that there is the local equilibrium value of the lattice constant. Hence, 
the lattice mismatch can be incorporated in the models for bandstructure cal- 
culations by defining the strain associated with it as a mismatch between two 
material layers 



e m = (a - a(r))/a°, (1) 

and by accounting for it in the strain-displacement relationships. In (1), o(r) 
and a are lattice constants of two different material layers, respectively, while 
r is the position vector responsible for tracking the interface. This aspect of 
mechanical effect contributions has been actively incorporated into the the- 
ory of bandstructure calculations since early 1970s, starting from fundamental 
works by Pikus, Bir, Rasba, Sheka and many others [15]. Secondly, semicon- 
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ductors are piezoelectric materials and the piezoelectric effect contributions to 
the overall properties cannot be ignored in bandstructure engineering, neither 
for hexagonal (wurtzite) structures (often due to the principle strain com- 
ponents) nor for cubic (zinc-blende) structures (often due to the shear strain 
components). Piezoelectrics represent anisotropic media and wave interactions 
in such media have been a topic of immense practical importance. While purely 
elastic waves in structures have been studied intensively for many decades, the 
study of coupled electromechanical interactions in anisotropic materials is of 
more recent origin. One reason for that lies with the fact that dealing with 
problems of coupled electroelasticity usually requires the development and 
implementation of effective numerical techniques [16]. Hence, the experience 
that has been accumulated in solving problems of coupled electroelasticity for 
piezoelectric structures becomes invaluable in the area of bandstructure engi- 
neering and modelling semiconductor quantum structures. As in the other ar- 
eas where modelling piezoelectric solids is an essential component, we consider 
the following general model (e.g., [16]), describing coupled electromechanical 
interactions in the Cartesian system of coordinates x = (xi, x 2 , xs) T : 

^ =V ' <T + F ' divD = G, E = -Vy?, (2) 

where u = (ui, u 2 , u 3 ) T is the displacement vector, cr = (cr^) is the stress, p 
is the density of the piezoelectric material, F and G are body and electric 
forces on the piezoelectric, if any, E and D are the electric field and electric 
displacement, and <p is the electrostatic potential. The conventional procedure 
applied in modelling LDSNs is based on the minimization of the purely elastic 
functionals (e.g., [4]), rather than on the solution of the fully coupled problem. 
The effect of coupling has been analyzed rigorously in a general setting in 
[18,19,16] (see also references therein), while in the context of nanostructure 
modelling it has recently been demonstrated that such an effect could be 
quite substantial [6,7]. The core component of our model for analyzing the 
properties of quantum dot nanostructures will be the equilibrium equations of 
the coupled theory of electroelasticity which are simplified in this case to 

daij/dxj = 0, divD = 0, (3) 

where the coordinate subindeces in the Timoshenko-Karman notions are ob- 
tained by changing 1 — > x, 2 — > y, z — > 3 in the tensorial representation above. 
The type of coupling between the mechanical and electric fields is determined 
by the type of the crystallographic symmetry of the material. In particular, 
for the WZ semiconductors we have 

&xx = Cne xx + Ci2£ yy + Ci 3 E zz — 613.^, <J xy = (c n — Cl 2 )e xy /2, 
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Cyy = ci 2 e xx + cnEyy + c\ 3 e zz — e 3 iE z , o yz = c^e yz — e^Ey, 

&zz = C\3\ £ xx + e yy) + C 33 E ZZ — G 33 E Z , <J ZX = C^£ zx — ei 5 E x , 

D x — £i5£zx + £i\E Xl D y = eiz,£y Z + €\\E yi 

D z = e 31 (e xx + e yy ) + e 33 £ zz + e 33 E z + P sp , (4) 

where and ta are piezoelectric and dielectric coefficients; P sp is the spon- 
taneous polarization. While for the WZ materials the built-in spontaneous 
polarization and the principal components of strain are main contributors 
to the piezoelectric effect contributions, in ZB materials it is the shear strain 
components that may contribute noticeably to the overall properties. For such 
materials we have the following constitutive relationships that couple (3): 

&xx C\\E XX -\- CyiEyy ~\~ C\2^xxi ®yy ^Yl^-xx CnEyy -\- CyiE ZZl 
®zz = C\2E XX + C\2E yy + C\\£ Z zi &yz = ^C^E yz — G\aE X) 
O "zx = ^44&zx Gi4-£7y, <7 X y = AC/nE X y C\^E Z i 

D x = euEyz + €uE x , Dy = eu£ zx + ^Ey, D z = e\±s xy + e 33 E z . (5) 

The issue of coupling via boundary conditions remains largely untouched in 
the area of modelling piezoelectric semiconductor nanostructures, in particu- 
lar when the wetting layer is taken into account. The formulation of correct 
boundary conditions in the latter case was discussed in [8]. Recall that the 
model we presented in [8] accounted for electron states with arbitrary kinetic 
energies in the wetting layer. Effectively, the model we derived for the quantum 
dot structure with wetting layer allowed us to demonstrate several important 
observations. In particular, electron states that correspond to a single quan- 
tum dot structure with wetting layer will asymptotically approach one of the 
two limiting situations: either "pure" quantum well states far away from the 
quantum dot region or zero in the case of a "pure" quantum dot state. This 
observation must be used, as explained in [8], for the formulation of general 
boundary conditions for the combined quantum-dot/wetting-layer structure. 

The resulting problem we deal with in this paper is a boundary value problem 
that is solved with respect to (u\,U2,u 3 , ip). From a mechanics point of view, 
the model is derived from a variational principle applied to the total potential 
energy which includes both deformational energy and piezoelectric field func- 
tional as described in [19,16]. Variational difference schemes developed in [16], 
as well as the finite element formulation developed for computations in this 
paper, follow from such a variational representation. Finite element method- 
ologies have been previously applied to bandstructure analysis in [3,17,20]. 
However, in these papers the contribution of piezolectric effect was not ac- 
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counted for. All works in this area we are aware of are based so far on the 
linear theory of elasticity. 

Taking into account piezoelectric effect contributions, in the subsequent sec- 
tions we will compare the results for bandstructure calculations and the predic- 
tion of optoelectromechanical properties of nanostructures that are obtained 
with linear and nonlinear strain models. 



4 Bridging the scales to the quantum effect level and exploiting 
the analogy with coupled models of structural mechanics 

Already today, quantum effects play an important role in many optoelec- 
tronic devices and structures. This trend will persist into the future as device 
miniaturization continues. Therefore, ideally the full bandstructure transport 
description is required. However, at present it is not possible in practice, in 
particular at the device/structure level. The reason is simple: transport should 
be computed with a many-particle Hamiltonian for the carriers and the atomic 
structures of the device/structure material. This is a task of enormous com- 
putational complexity, not feasible to complete today. Hence, some simplifica- 
tions need to be made. 

One approach is to account for quantum mechanical (and statistical) effects via 
the Wigner-Boltzmann model. Although such a model also involves substantial 
computational difficulties, it allows us to construct a hierarchy of the macro- 
scopic models in a way similar to those involving fluid dynamics problems 
and semiconductor device theory where the continuity (fluid-like) analogy is 
used for the model classification (e.g., [21] and references therein). At present, 
most of the quantum corrected macroscopic models are in their infancy as they 
are usually not able to adequately include interactions between the electrons 
and other particles. New efforts in this direction are currently being under- 
taken by a number of authors (e.g., [22]). Since the problem we are addressing 
is a multiscale problem, a natural way to approach its solution in the above 
framework could be to apply a domain decomposition technique. For example, 
one can use the quantum mechanical approach in the regions where quantum 
mechanical effects are dominant and use continuum-like (e.g., hydrodynamic) 
models in other regions. However, the issue of coupling such models, e.g., via 
an interface condition of a typical domain decomposition methodology or by 
using other techniques, is far from trivial and remains largely open. 

In this paper we follow another route. While the application of ab initio and 
atomistic methodologies are inheritably problematic from a computational 
complexity point of view, we resort to averaging procedures over atomic scales. 
This can be achieved by a variety of procedures, including various empirical 
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tight-binding, pseudopotential, and k • p approximations. In what follows, we 
focus on the latter approximation as a tool for averaging over atomic scales. 
The procedure stems from the original work by Luttinger-Kohn and is based 
on the effective mass approximation and the subsequent development of the 
k • p theory. As with any model, the one we develop here relies on a set of 
assumptions some of which may not be always fulfilled. For example, a typical 
assumption of the kp theory that potentials change slowly on the length scale 
of the lattice constant could be questionable for Metal-Oxide-Semiconductor 
Field-Effect (MOSFE) devices, e.g. Nevertheless, in bandstructure calculations 
of LDSNs, the theory provides a remarkably flexible tool. Furthermore, while 
we do not address this issue in detail in this paper, it is worthwhile mention- 
ing that a recent refined approximation of the Luttinger-Kohn Hamiltonian, 
known as the Burt-Foreman correction, has been developed and tested (see 
examples and further details in [23]). This correction allows us to put the effec- 
tive mass theory formalism related to the behaviour of the envelope functions 
across interfaces on a much more rigorous mathematical foundation. 

For the benefit of the reader, we recall the main premises of the k-p approxima- 
tion. Although a number of methodologies quoted above (such as tight-binding 
and pseudo-potential) can provide us with the global dispersion relationships 
over the entire Brillouin zone for the bulk material, to know the main elec- 
tronic characteristics of the semiconductor, such as wave functions, we need 
only the dispersion relationship over a small wave vector k around the band 
extrema [24]. Indeed, it is well known that most processes in semiconductors 
take place near the top of the valence band and at the bottom of the conduc- 
tion band [25]. The wave functions, characterizing wave propagation in the 
box volume Q = {(x,y,z) : < x < L x ,0 < y < L y ,0 < z < L z ,}, are 
assumed to be in the form 

$W(r) = (L x L y L z )~ 1/2 exp[i(k x x + k y y + k z z)\ = 

(L x L y L g )-^ 2 exp(ik-r), (6) 
where the allowed values of k 



k= (2irl/L x ,2irm/L y ,2im/L z ), I, m, n = 0, ±1, ±2, ... (7) 

form a 3D k space with the origin denoted as T point. Based on this as- 
sumption, the application of Bloch's theorem in the context of semiconductor 
crystals leads to the following (Bloch) representation of the wave function in 
a crystal 

$ nk (r) = u nk exp(ik • r) (8) 
with -u nk being periodic. Since it is easier to find approximate solutions for the 
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function w ra k (assumed to be slowly varying over a small region of k space) 
than for $ n k [25], we attempt to write the Schrodinger equation in terms of 
w n k. In doing so, the derivatives of the momentum operator p = —ihV acting 
on the plane waves are simplified to Kk, leading to a simplification of the 
Schrodinger equation for the Bloch functions of the crystal, $ n k, where the 
terms depending on k are treated as perturbations away from the solution at 
k = 0. The presence of the operator kp in the resulting expression renders the 
name of the underlying local methodology. In this way, the standard parabolic 
approximations of the bands are refined locally in the k • p theory. This is 
important for the conduction band, but even more so for the valence band 
due to degenerating light and heavy holes at the T point. 

It is worthwhile mentioning that another major reason for this particular 
choice of averaging lies with the fact that the k ■ p treats the bandstructure 
in a continuum-like manner, albeit allowing to incorporate many important 
effects acting at different scales which is very important for the problems like 
ours. This includes mechanical and electromechanical effects. As strain and 
piezoeffect are key contributors to changes in optoelectromechanical proper- 
ties, the models described in the previous section must be incorporated in 
the bandstructure calculation. In order to do that, we first remind the reader 
that the accuracy of approximations based on the k • p theory depends on the 
functional space where the envelope function is considered. In fact, we put 
subbands within conduction and valence bands of the semiconductor material 
into correspondence to the basis functions that span such a space. The num- 
ber of such functions varies in applications, but typically ranges from 1 to 8 
to ensure computational feasibility of the problem. For WZ materials, 8 func- 
tions should usually be included due to spin-orbit, crystal-field splitting, as 
well as conduction/valence band mixing effects. This corresponds to 6 valence 
subbands and 2 conduction subbands that account for spin up and spin down 
situations. Given that the Hamiltonian in the k • p theory can be represented 
as 



the problem at hand can be formulated as an eigenvalue partial differential 
equation problem 



The standard Kohn-Luttinger representation of the Hamiltonian in the form 
of (9) is well documented in the literature and we refer the reader interested 
in details to books by G. Bastard and J. Singh, e.g. [24,26]. The problem (10) 
should be solved with respect to eigenpair (\l>, E), where E is the electron/hole 




(9) 



= E^. 



(10) 
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energy and \1/ is the wave vector with dimensionality of the functional space 
chosen. For example, in the situation discussed above we have 

* = (4, 4 , Vi 4, 4, 4) T , (ii) 

where the subindex S denotes the wave function component of the conduction 
band and ipx = (\X > | f) is the wave function component that corresponds 
to the X Bloch function of the valence band when the spin function of the 
missing electron is up. In the Hamiltonian representation (9) 7i is the energy 
functional, defined either by the standard Kohn-Luttinger Hamiltonian (as 
mentioned above) or by the Burt-Foreman modification [23]. It represents the 
kinetic energy plus a nonuniform potential field V and other effects contribut- 
ing to the total potential energy of the system, as specified further in Section 
4. Other notations are standard: h is the Planck constant, m is the free elec- 
tron mass, r = (x\ = x, X2 = y,x^ = z), while the superindeces (m,n) are 
used to denote the basis of the space for the wave function, which in the case 
discussed above (6 valance and 2 conduction subbands) would lead to an 8X8 
Hamiltonian. 

Four equations (3) of coupled electromechanics and the eigenvalue PDE prob- 
lem (10) constitute a coupled Schrodinger-Poisson model that provides a gen- 
eral framework for incorporating nonlinear effects. Similar to the top-down 
approach, where coupling procedures are well established, in the bottom-up 
approach applied here the coupling between Schrodinger and Poisson equa- 
tions is essential, in particular at the level of quantum device modelling. 



5 Coupling strain with bandstructure calculations in the varia- 
tional formulation 

Since the fundamental work [15], strain effects firmly took their place in the 
models for bandstructure calculations of semiconductor structures. Many im- 
portant examples of such calculations based on strained Rasba-Sheka-Pikus 
Hamiltonian have recently been provided in the context of low-dimensional 
semiconductor nanostructures, including quantum dots (e.g., [4]). However, 
all current models we are aware of have been based on linear theories. The 
first question to ask is whether material nonlinearities, expressed by stress- 
strain relationships, may become important for such calculations. The usual 
argument for using the linear stress-strain relationships in this field of ap- 
plications is based on the fact that strain is indeed of order of magnitudes 
smaller of the elastic limits. Although this argument may fail at the device 
level modelling, it remains valid for the structures of interest in the present 
paper. At a more fundamental level, however, elastic and dielectric coefli- 
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cients may become nonlinear due to geometry of the structure which will lead 
to nonlinear stress-strain relationships. To address this issue rigorously, one 
has to look critically at the standard Keating model where the parameters 
of the model are unit-cell dimension dependent. In the general case, we may 
need to account for many body interactions of higher order to reflect asym- 
metry of the interatomic potential. This issue is outside of the scope of the 
present paper. Instead, we will focus on a consequence of this issue that leads 
to geometric nonlinearities. Indeed, one of the major drawbacks of the cur- 
rent models for bandstructure calculations is that they are not able to resolve 
adequately strain nonhomogeneities due to the application of the original rep- 
resentation of [15] based the infinitesimal theory with Cauchy relationships 
between strain and displacements. As we demonstrate in the next section, 
this approximation is inadequate when we have to deal with geometric irreg- 
ularities of low-dimensional semiconductor structures. In the latter case, by 
using the variation of deformation be^ which is induced by the variation in 
displacements Sui, we describe the new position £j = xi + Ui of the material 
particle (with initial coordinates Xi) after deformation. This leads naturally to 
the formulation of the problem where the general nonlinear Green-Lagrange 
relationship for strain can be accounted for. As we pointed out, the varia- 
tional formulation of the problem at hand has been used before in the context 
of finite element implementations (e.g., [3,17,20]). Two new features that have 
been developed in this paper include (a) consistent coupled treatment of the 
piezoelectric effect and (b) the ability to incorporate geometric nonlinearities 
into the model. As we demonstrate in the next section, the influence of these 
features on the bandstructure calculations, and therefore optoelectromechan- 
ical properties of the modelled structures, could be substantial. 

The case of weakly coupled piezoelectricity was analyzed previously in [18] 
where convergence results in Sobolev classes of generalized solutions were es- 
tablished. Indeed, in some special cases, in particular for simple boundary 
conditions, the problem can be addressed semi-analytically. In the context of 
quantum dots, this has been done in [6] by using Green's function approach. 
Strongly coupled problems of piezoelectricity have been analyzed previously 
from variational scheme perspectives in [16]. In the latter case, the steady- 
state formulation given by (3) should be understood in a variational sense 
as: 



where both piezoelectric stress and the nonlinear part of strain are taken into 
account via constitutive relationships with the total variation of deformation 
given by its linear and nonlinear parts 5e = 5el + Se^. 

The above problem is coupled to the eigenvalue PDE problem, also understood 




(12) 



v 
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in a weak sense. That is we seek the solution to the following problem 



j_2 

-> min, = / (W) T H {m ' n) Wdv - E [ ^ T qdv, (13) 



t J 



where the Hamiltonian is given by (9). As usual in the k • p theory, formal 
representation can be reduced to the sum of constant and k-dependent energies 
(H = H + V): 



H = H + j2H t , (14) 

i=l 

where H , derived from the standard Kane Hamiltonian at k = (e.g., [24]), 
accounts for the spin-splitting effects; Hi is the contribution due to the kinetic 
part of the microscopic Hamiltonian unit cell averaged by the respective Bloch 
functions, S, X, Y, or Z; H 2 is the strain-dependent part of the Hamiltonian, 
and H 3 is the energy of unstrained conduct ion/ valence band edges. 

In what follows we consider an example of modelling nanostructures based 
on the model described above. The main emphasis is given to a pyramidal 
quantum dot residing on a wetting layer. In [8] it was demonstrated that the 
electronic states in the wetting layer may influence electronic states in the dot 
and vice versa and, hence, the thin layer on which the quantum dot resides 
cannot be excluded from the computational domain as it is often done in the 
literature. This brings additional difficulties in the computational implemen- 
tation of this multiscale problem. 



6 Computational experiments on quantum dot nanostructures 



Our representative example in this section concerns InAs/GaAs quantum dot 
structures of pyramidal shape. Such self-assembled structures are grown exper- 
imentally via Stranski-Krastanov methodology and have been studied before 
with simplified models of linear elasticity (e.g., [1,5]). As we have already men- 
tioned, the present study has a number of new features. Firstly, we take the 
full electromechanical coupling as well as the wetting layer into account. Sec- 
ondly, we analyze the influence of nonlinear strain components in modelling 
such quantum dot structures. 

All our computational examples below are for the quantum dot structure pre- 
sented in Fig. 1 where geometric dimensions are also given. Pyramidal shapes 
of quantum dot structures, such as the one we analyze, have been confirmed by 
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Wetting layer 

Fig. 1. Geometric dimensions of the representative pyramidal quantum dot. 
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Fig. 2. Energy levels of the conduction and valence bands under the linear strain 
approximation . 

experimental techniques, including high resolution electron microscopy. The 
entire structure, consisting of the InAs dot sitting on the wetting layer, is 
embedded in a (spherical) GaAs matrix. In such semiconductor materials one 
expects changes in the optoelectromechanical properties due to strain effects. 
However, up until now, such changes have been quantified with linear theo- 
ries only. Hence, as the first step, in Fig. 2 we present energy levels for both 
conduction and valence bands by using the conventional methodology based 
on the linear approximation. This result is given at the center of the dot 
(x — 0, y — 0) along vertical z-axis. 

A straightforward generalization of the linear theory in our context is to ac- 
count for the large deformation gradient Vu^. This gradient is responsible for 
the dominant nonlinear strain effect due to the lattice mismatch in the growth 
direction. While in the x and y directions the strain components will be iden- 
tical in this case to the von Karman type model, in the z-direction they differ. 
As all the shear strain components are assumed in this case [Case I] to be 
linear, the Green-Lagrange strain components in this model take the form: 

_ du ± | If du 3 \ 2 _ du 2 , 1 f du 3 \ 2 

£xx ~ dx + 2\dx J ' £yy ~ dy + 2\dy ) ' 
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Fig. 3. Energy levels of the conduction and valence bands accounting for the large 
deformation gradient in the growth direction. 
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Fig. 4. Quantifying nonlinear contributions in the band edge potentials, accounting 
the large deformation gradient in the growth direction. 



du 3 1 / du 3 
~dz + 2 \!h 



(15) 



In this case, closer to the base of the dot, the behaviour of the energy levels 
are the same as in the linear case. However, the situation changes at the tip 
of the dot as can be seen from Fig. 3. As before, the result is presented at the 
center of the dot (x — 0, y — 0) along vertical z-axis. 

Next, we have analyzed and quantified the difference between this nonlinear 
case and the linear case, conventionally used in these calculations. In Fig. 4 we 
present the difference in the band edge potential for the conduction and valence 
bands. This result demonstrates that nonlinear strain contributions could be 
substantial in calculating energy levels, and hence in predicting optoelectronic 
properties of nanostructures. In our case these contributions are particular 
pronounced for the conduction band. 

The comparisons indicate the maximum deviation from the linear theory near 
the top vertex of the pyramid of 60nm high, not at the wetting layer. In- 
deed, due to the lattice mismatch at the wetting layer and pyramid interface, 
components like (dui/dx) 2 and (du 2 /dy) 2 will become important, rather than 
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Fig. 5. Biaxial and isotropic strains at y=0. 




Fig. 6. Biaxial and isotropic strains at z=1.5nm 

(du 3 /dx) 2 or (du 3 /dy) 2 which are due mainly to the Poisson effect in this 
problem. 

These differences between the results produced with the linear and nonlinear 
models are hard to quantify without the band edge potentials, as presented 
above. Indeed, we have calculated a number of other characteristics, including 
biaxial strain 

£b = (£ll - E-a) % + (^22 - £33) 2 + (^33 - ^llf (16) 

and isotropic strain 

Si = s n + £22 + ^33- (17) 

In Fig. 5 we present biaxial and isotropic strains in the y=0 plane. 

These characteristics are shown also just above the wetting layer, at z=1.5nm, 
in Fig. 6. In both of the above situations we chose to present the results of 
calculations produced with the conventional methodology, while noting that 
calculations obtained with nonlinear contributions look very similar. 

In Fig. 7 we present the distribution of the piezoelectric potential in the quan- 
tum dot under consideration (as well as its projection from the top). Recall 
that results with uncoupled or semi-coupled models usually demonstrate max- 
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Fig. 7. Distribution of the piezoelectric potential in the quantum dot. 



Fig. 8. Ground state of the quantum dot nanostructure. 

ima of the piezopotential outside of ZB quantum dot structures only, in partic- 
ular when the dot is truncated. This is, of course, not the case for the coupled 
model applied in the current situation. The local extrema of the piezoelectric 
potential near the QD top can be well reproduced, as demonstrated by Fig. 
7. As expected, the distribution is symmetric at the bottom and at the top of 
the pyramid. We note also that predicting optoelectromechanical properties of 
hexagonal (WZ) materials usually requires to account for an additional effect 
of spontaneous polarization, typically negligible in ZB materials. As described 
in Section 2, our model is capable of dealing with both types of the materials. 

Next, we have calculated eigenstates of the structure. The ground state, pre- 
sented in Fig. 8, has been calculated accounting for strain and piezoelectric 
effects. As seen from the figure, the state is fully confined. 

In Fig. 9 the next four eigenstates of the nanostructure are presented. As 
confinement visualization looks similar for both linear and nonlinear cases, we 
shall provide calculated numerical values that demonstrate the influence of 
piezoelectric effect as well as nonlinear contributions on electronic states. 

Before doing that, we recall that the only nonlinear model we considered so 
far was the model accounting for the large deformation gradient in the growth 
direction. We have also analyzed the case where we account for the large 
deformation gradients V«i and V«2, responsible for the dominant nonlinear 
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Fig. 9. Electronic confinement in the quantum dot nanostructure: the 2nd, 3rd, 4th, 
and 5th eigenstates. 

strain effects due to the lattice mismatch in the x — y plane, normal to the 
growth direction. In this case [Case II] the Green-Lagrange strain components 
are: 




_du 2 if duA 2 1 / du2 \ 2 
£yy ~ dy + 2\dy + 2\dy 



dus 1 ( du\\ 2 1 / du 2 x 



- 8: + 2[l>t) + 2[W) ' (18) 

while all the shear strain components remain linear. As expected in this case, 
the results show nonlinear effect contributions due to the nonlinear relaxation 
of the interfacial strain between the InAs wetting layer and the GaAs matrix. 

Next, we implemented the full nonlinear strain model [Case III] in the devel- 
oped finite element code, based on the general Green-Lagrange strain 



e = I(F T F-I), (19) 

where F is the deformation gradient and I is the identity matrix. In Figs. 10- 
12 the results are presented along the z-axis (y=0) for the three cases: x=0, 
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Fig. 10. The band edge potential with the full nonlinear strain model (x=0, y=0). 
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Fig. 11. The band edge potential with the full nonlinear strain model (x=30, y=0). 

30, and 60 nm, respectively. Although differences with calculations obtained 
with the linear model are clearly observed, we quantify them in Table 1. The 
first column in the table is the number of the corresponding eigenstate. The 
second column gives eigenstates obtained with the linear theory, accounting 
for strain only; the third column gives eigenstates obtained with the linear 
theory, accounting for strain and piezo effects. The remaining columns give 
the results obtained with the nonlinear models based on Cases I, II, and III, 
respectively, as described above. All values are given in eV. 

In all finite element computations reported here we applied tetrahedral ele- 
ments with quadratic Lagrangian interpolation function. The global conver- 
gence was analyzed by refining the finite element mesh and ensuring that the 
difference between the last two subsequent refinements is negligible in L 2 . For 
the solution of the discretized equations we used GMRES with incomplete 
Cholesky factorization. The maximum number of iteration was set to 500 and 
the L 2 error tolerance was set as 1 x 10~ 6 . The pyramid was embedded into a 
sphere. The loading conditions were simulated by subjecting the outer surface 
of the sphere to Dirichlet boundary conditions in displacement. The pyramid 
structure was allowed to equilibrate under the combined effect of lattice misfit 
induced strain and piezoelectric polarization. 

Finally, we calculate the bandstructure with our finite element method code 
and compare the band edge potential due to the above three nonlinear strain 
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Fig. 12. The band edge potential with the full nonlinear strain model (x=60, y=0). 



Eigenstate # 


Lin/strain 


Lin / str ain+piezo 


Case I 


Case II 


Case III 


1 


0.7122 


0.6968 


0.6811 


0.6913 


0.6767 


2 


0.8345 


0.8084 


0.8008 


0.8070 


0.7994 


3 


0.8404 


0.8115 


0.8046 


0.8095 


0.8027 


4 


0.8511 


0.8272 


0.8248 


0.8267 


0.8244 


5 


0.8658 


0.8350 


0.8321 


0.8345 


0.8316 



Table 1 

The influence of strain, piezoeffects, and nonlinear contributions on eigenstates of 
the structure. 




(nm) 



Fig. 13. Quantifying the difference between linear and nonlinear models in the pres- 
ence of piezoeffect (x=0, y=0). 



models in the presence of piezoelectricity. The difference for each of these three 
cases of band edge potentials with respect to that of the linear case is computed 
and plotted along the z-axis for (a) x — 0, y — 0, (b) x = 30nm, y — 0, and (c) 
x = 60nm, y = in Figs. 13-15, respectively. It is clear from the plots that 
the fully nonlinear strain model (Case III) captures the large inhomogeneous 
strain-induced band edge potential as a cumulative effect of the potentials for 
Cases I and II. The band edge potentials along the presented [001] direction (z- 
axis) for each of these three cases have been computed with the fully coupled 
model. 
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Fig. 14. Quantifying the difference between linear and nonlinear models in the pres- 
ence of piezoeffect (x=30, y=0). 
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Fig. 15. Quantifying the difference between linear and nonlinear models in the pres- 
ence of piezoeffect (x=60, y=0). 



7 Conclusions 



In this contribution, we analyzed two fundamental approaches in bridging 
the scales in mathematical models for the description of opto electromechani- 
cal properties of nanostructures. We demonstrated that the conventional ap- 
plication of linear models to the analysis of properties of nanostructures in 
bandstructure engineering could be inadequate. By accounting consistently 
for the piezoelectric effect, we considered three nonlinear strain models and 
analyzed contributions of quadratic nonlinear terms induced by the deforma- 
tion gradient in the growth direction, as well as in the plane normal to that 
direction. The core of our model was based on the Shrodinger-Poisson sys- 
tem which was presented in the variational form and implemented with finite 
element methodology. In this framework, we presented also the full nonlin- 
ear strain model and quantify the differences between the conventional linear 
and developed here nonlinear models for bandstructure calculations. Gener- 
alizations of the existing models and the examples provided for quantum dot 
nanostructures emphasized the importance of coupled effects in predicting 
optoelectromechanical properties of such structures. 
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